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The conformally flat families of initial data typically used in numerical relativity to represent 
boosted black holes are not those of a boosted slice of the Schwarzschild spacetime. If such data are 
used for each black hole in a collision, the emitted radiation will be partially due to the "relaxation" 
of the individual holes to "boosted Schwarzschild" form. We attempt to compute this radiation by 
treating the geometry for a single boosted conformally flat hole as a perturbation of a Schwarzschild 
black hole, which requires the use of second order perturbation theory. In this we attempt to mimic 
a previous calculation we did for the conformally flat initial data for spinning holes. We find that the 
boosted black hole case presents additional subtleties, and although one can evolve perturbatively 
and compute radiated energies, it is much less clear than in the spinning case how useful for the 
study of collisions are the radiation estimates for the "spurious energy" in each hole. In addition 
to this we draw some lessons on which frame of reference appears as more favorable for computing 
black hole collisions in the close limit approximation. 
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There is significant interest in obtaining waveforms for the gravitational radiation produced in the collision of black 
holes. Progress is being made on this problem both using supercomputers |l|] and perturbative calculations ||. One of 
the open issues is what families of initial data are appropriate to represent the collision of two black holes, especially 
when the latter are not far away from each other. The state of the art of numerical simulations suggests that for 
some time we may not be able to start simulations with the black holes at a sufficiently large separation, such that 
^ ! one can assume a simple linear superposition will work. This leaves open the issue of how much "spurious radiation" 
I , is one introducing in the various proposals for superpositions in the non-linear regime. Bowen and York || (and with 
J~ I ■ a different set of boundary conditions more recently Brandt and Briigmann [0 ) studied the problem of giving initial 
data for boosted and spinning holes in such a way that a superposition is possible. They assume the spatial metric 
is conformally flat and as a consequence one can superpose the extrinsic curvatures for the two holes and still solve 
' the momentum constraint. One then proceeds to find a conformally flat spatial metric for the superposed holes by 
solving a nonlinear elliptic partial differential equation. The procedure achieves superposition at the price of assuming 
conformal flatness of the three metrics, which is not generically possible, and more importantly, is not possible in 
situations of interest. For instance a single spinning black hole, described by the Kerr solution, is not known to admit 
conformally flat spatial slices pi. 

A similar situation develops for the case of a single boosted black hole. The initial data constructed by Bowen and 
York or Brandt and Briigmann do not correspond to those one would find on a boosted slice of the Schwarzschild 
spacetime. In this paper we will refer to these families as "conformal boosted black hole" Therefore if one evolves 
these families of data one should find that the black hole "settles down" to a Schwarzschild form through the emission 
of gravitational radiation. The original purpose of this paper was to study the emitted radiation by treating the 
conformal boosted hole as a perturbation of a Schwarzschild black hole. In a previous paper we had carried out a 
similar discussion for the case of a single spinning conformally flat hole H . We will see that the boosted case is much 
more subtle than anticipated. We will be able to evolve the spacetime, but questions will remain about the usefulness 
of the results obtained, at least for the original purpose of gaining insight into the spurious radiation content of data 
for black hole collisions for interesting ranges of parameter values. 

The organization of this paper is as follows. In the next section we will discuss perturbations of a boosted black 
hole. We recall that the £ = 1 even modes are pure gauge and therefore all the physics of interest takes place in the 
£ = 2 modes, which can be treated easily up to second order in perturbation theory. In section III we will discuss 
the perturbative evolution and the amount of radiation produced. We end with a discussion of the results and their 
implication for the choice of frame of reference one uses in perturbative evolutions of black hole collisions with linear 
momentum. 
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II. CONFORMAL BOOSTED BLACK HOLE AS A PERTURBATION OF SCHWARZSCHILD 



A. Initial data in the conformal approach 

The families of initial data that we will consider in this paper are obtained via the "conformal approach" to the 
initial value problem in general relativity. In it, one assumes the metric to be conformally flat g a b = 4- 'Sab and defines 
the conformal extrinsic curvature K a b — <f> 2 K a b- In terms of these variables the initial value constraint equations 
(assuming maximal slicing TrK — 0) read, 

V a K ab = (1) 
1 K ab K >, 

where all the derivatives are with respect to flat space. 

One can construct || solutions to the first set of equations (momentum constraint) for a single black hole centered 
at R — 0, with linear momentum P a , 

■~ 3 

K ab = [ 2P (anb) - {Sab - n a n b )P c n c ] . (3) 

where R is a spherical radial coordinate and n a a radial unit vector, and both are defined in the fiducial flat space 
that one obtains setting the conformal factor to unity. 

Without loss of generality we may assume that P a points along the positive z-axis, and has magnitude P. If we 
write K a b in spherical coordinates, the only non vanishing components are 

-» 3P 
K RR = -^j cos(0) 

^ 3P ~ 

K m = -^ sin ( 9 ) = k or 

— 3P 

Kee = — — cos(0) 

K u = sm{6) 2 K ee (4) 

If we write the extrinsic curvature in terms of tensor spherical harmonics (^] , we see that it consists of a pure £ = 1 
even term. This is reasonable, since the presence of momentum in an initial slice configuration is determined by the 
presence of a "dipole" term asymptotically, that makes the ADM integral 

p i = ^zl K ijd 2 S j (5) 



nonvanishing. 

From (0), we have, 
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KabK =^l( 1 + 2cos20 ) ( 6 ) 



We now need to solve (^), and this involves imposing boundary conditions on <f>. One possibility is that given by 
Bowen and York ||[]. In this case one chooses a certain constant a, solves (||) for R> a, and imposes 

ff + ^ = 0for*=a, (7) 

and 



1 Alternatively, we may use a Brill-Lindquist type boundary condition j?| and as a result one obtains the "puncture" solutions 
of@. 
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> , lim = 1. 



(8) 



It should be clear that, in addition to the choices we originally made, there are other choices made in constructing 
the solutions. As suggested by Bowen and York, a possibility that can be considered appealing, especially when 
one has multiple black holes, consists in requesting that both the conformal factor and the extrinsic curvature be 
symmetrized through the throats of the holes. This was implemented explicitly by Cook Q in his numerical work 
on constructing solutions for the conformal factor for multiple black holes. The symmctrization procedure yields 
non-trivial results even for a single black hole as we are considering here |lC| . We have chosen, for simplicity, not to 
symmetrize the extrinsic curvature. However, the boundary condition (Q) automatically ensures that the conformal 
factor be symmetric (provided one chose a symmetric extrinsic curvature). One could perform other choices and the 
problem is simple enough that — at the level of approximation we are working — it can always be solved. Perhaps 
the choice we make here is not the most aesthetically appealing to some readers, but experience has shown that 
symmetri zing or not symmetrizing does not change significantly the amount of radiated energy in head-on black hole 
collisions . Therefore our choice should not crucially influence the central conclusions we are attempting to obtain 
about the conformally flat black hole solutions. It should be noticed that with our choice of non-symmetric extrinsic 
curvature, condition (Q) implies that the surface where it is applied is an extremal surface (not necessarily minimal). 

Even for the simple form of (Q), in general we cannot solve (|^) exactly, and one has to resort to numerical methods. 
In the present analysis, however, we will be interested in solutions for "small" P. Since for P = the solution is 

(O) =1 + | (9) 

we may solve (^|) iteratively replacing 

= 0(0) (R) + P 2 ^ (R, 9) + P 4 0( 4 ) (R } e) + ■ ■ ■ (10) 



in (J2J) , and imposing 



°" + ^-(f> il) =0 for R = a, lim = 0. for i ^ 0, (11) 



as 



The form of @ further suggests that we expand <p^ in Legendre polynomials Pi(cos9), so that 0( 2 ) may be written 

<t> = ^\R) + P 2 [4 2) (i?)P o (cos0) + 4 2) (^)^(cos0)] +0{P i ). (12) 
Solving for the coefficients taking into account the boundary condition (j?]) we get, 
W fm _ (R + a) 5 ~R 5 , 13 



32a 2 (R + a) 5 512aR 



(2 ) fm _ 75R 6 + 291aR 5 - 650a 2 i? 4 - 3400a 3 i? 3 - 4800a 4 i? 2 - 2925a 5 i? - 669a 6 

400i? 3 (i? + a) 5 



21a r , , , / „ 121a 

+ ? n(2a) - \n(R + a) } + (13) 

20i? 3L y ' 128i? 3 V ' 

Notice that (jxp {R) falls off only as R~ 2 for large R. 

For completeness, we also present the solution one obtains if one chooses the "puncture" boundary condition 
considered by Brandt and Briigmann B recently, 

,(2l,™ Af 4 + 10 M 3 R + A0M 2 R 2 + 80 MR 3 + 80 i? 4 . „ 

^ {R) = 8 M(2R + Mf (M) 

, ( 2i,™ 120 R 5 + 768 R 4 M + 1078 R 3 M 2 + 658 R 2 M 3 + 189 RM 4 + 21 M 5 . , 

fflo '(R) = = (15) 

20i? 2 (2i? + Af) 5 

(ln(2i? + M)-ln(M))M 

" 40i^ 

but we will only consider the solution (O) in the rest of this paper. 
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To obtain an initial data set appropriate for a perturbative evolution we proceed in two steps. If, in an ADM type 
decomposition, we choose our shift functions iVj = 0, we have, 

ds 2 = g tj dx l dx 3 - N 2 (dt) 2 

= -2NK i:j (16) 



dt ~ 



where, for the initial slice, that we may take as t = 0, the 3-metric gtj and extrinsic curvature Kij are given by 
the above construction. We next change the conformal spherical r adial coor dinates R, to a "Schwarzschild" radial 
coordinates r, with R = (y/r + \Jr — 2M) 2 /4, and choose N = wl — 2M/r, so that we recover the Schwarzschild 
metric for P = 0. The perturbation treatment refers to this last form of the metric. It should be noticed that the 
extrinsic curvature to consider in (lm should be the physical one obtained dividing by the conformal factor squared 
the conformal extrinsic curvature (|). 

Carrying out all these steps we find that the initial data has the following multipolar components: to zeroth order 
in the momentum, we only have £ = components; to first order in the linear momentum, only £ = 1 contributions; 
to order 0(P 2 ), we have £ = 0, 2 contributions. All contributions are even-parity. 



B. Multipolar decomposition of the initial data : the i = 1 contributions 



The contribution to the three-metric of zeroth order in P is just the Schwarzschild solution. The first apparently 
non-trivial contribution is given at 0(P 1 ) and corresponds as we discussed in the previous subsection to an I = 1 
multipole. Let us analyze the I = 1 perturbations of a spherically symmetric spacetime. In order to do this we use 
the traditional Regge- Wheeler |q] decomposition. One starts with a background metric written as, 



Mdx^dx" = -(1 - 2M/r) dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 d9 2 + r 2 sin 2 9 dtp 2 



(17) 



For axisymmetric perturbations, the general 



h 
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I = 1 even parity terms take the form || 

(1 - 2M/r)H (t,r) cos (9 
Hi(t,r) cos 9 

(1 - 2M/ r y 1 H 2 (t,r) cosfl 

— ho(t, r) sva9 

— hi(t, r) sin6> 
r 2 K(t,r) cos6» 

r 2 sin 2 9K(t,r) cos 9 



(18) 



The first order metric perturbation coefficients are not uniquely defined, but may be changed by "gauge transfor- 
mations" of the form M 



h {i) 



9^ifip£(i) ff^pC^i)^ 9pv£,(\)i[i 



(19) 



where the gauge 4-vector £^ is arbitrary, except for the requirement of axisymmetry. 



In particular, the 



1 even parity coefficients transform as 

,aM (i,r) , 



H (t,r 
Hx(t,r 
H 2 (t,r 

h Q (t,r 
hi(t,r 
K(t,r 



H (t,r) 
Hi(t,r) 
H 2 (t,r) 
h Q (t,r) - 
hi(t,r) - 



2M 



dt r{r-2M) 
r -2M 8M (t,r) r 



M x {t,r) 

dMx{t,r) 



r 

2M 



r(r - 2M) 
r-2M 



dr 



r - 2M dt 
dMx{t,r) 



r 
2 



2M 



K(t,r) - -Mi(t,r 
r 



M (t,r) 
Mx{t,r) 

2M 2 (t,r 



dr 

2 dM 2 (t : r) 



dt 

2 dM 2 (t,r) 
dr 



(20) 
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where the functions Mi are arbitrary. 

One can use this gauge freedom to go to a restricted gauge in which ho = hi = K = 0. This gauge is not completely 
fixed. One still can perform gauge transformations of the form, 

where f(t) is arbitrary. 

An interesting result, is that in this gauge it is straightforward to find the general solution of the linearized Einstein 
equations for £ = 1 even parity perturbations. The result is, 

1 (r-2A/) 2 dt 

H 1 F | ^ d2p i 

3(r - 2M) 2 1 3M(r - 2M) 2 dt 2 

H ^V^y Fl{tl (22) 

where F\(t) is an arbitrary function. Remarkably, one can show that this solution is pure gauge. Choosing f(t) = 
—Fi(t)/(6M) in (|2l| ) leads to vanishing gauge transformed quantities. 

We therefore see that the £ = 1 perturbations are pure gauge. This result was first noticed by Zerilli |l2|. This was 
to be expected in physical grounds since one could always imagine setting coordinates boosted in such a way that the 
black hole would not move. It has the further implication that we can compute the radiated energy by computing 
the £ = 2 perturbations in a gauge in which the £ = 1 perturbations vanish, and studying their evolution. 



C. Multipolar decomposition of the initial data : the I = 2 contributions 

The relevant £ = 2 perturbations are of second order in our perturbation parameter, P. In principle, when one 
works out second order perturbations of a given metric, the evolution equations one gets have the general form of a 
linear operator (similar to the one that acts in first order giving rise to the Zerilli equation) acting on second order 
quantities, equal to a quadratic "source" term formed with the first order perturbations [13]. The source term is 
complicated and is delicate to handle numerically when evolving the perturbations. A place where this was explicitly 
done was for instance in the evolution of boosted black hole collisions to second order pj| . The calculations are 
lengthy and complicated. Fortunately, in our case one can proceed in a different way. We have just shown that 
there is a gauge in which the first order perturbations (the £ = 1 ones) vanish. Therefore in that gauge one can 
write a second order Zerilli equation that is source-free. Moreover, the linear portion of the second order perturbative 
equation is exactly the same as the first order perturbative equation (the Zerilli equation), and this equation can be 
written in terms of quantities that are gauge invariant. We notice, however, that eliminating the first order £ = 1 
terms through a first order gauge transformation introduces second-order changes in the metric that are not a pure 
second order gauge transformation, and must be taken into account. One must then be careful in making sure that 
the initial data for the £ = 2 perturbations used in the Zerilli equation corresponds precisely to that gauge. In other 
words, we need to carry out a first order gauge transformation on the initial data that provides a new initial data 
corresponding to a gauge where the £ = 1 perturbations vanish. Since all perturbations satisfy equations that are 
second order in time, this requires that the £ = 1 terms of the metric, and their first time derivative vanish on the 
initial slice. If we consider (|2(]), this requires that the gauge vector components Mi are such that both the left hand 
sides in (po|), and their first time derivatives, vanish when evaluated at t = 0. This, in principle requires only the 
knowledge of the metric functions on the right hand side of (pp|), and their first time derivative, at the fiducial time 
t — 0. It turns out, however, that we also need second order time derivatives, 

to implement the second order gauge transformation required to obtain the £ = 2 initial data. These second 
order time derivatives may be straighforwardly evaluated from the corresponding Einstein equations for the £ = 1 
perturbations. 

For the particular case in question, the £ = 1 initial data (and first time derivative) is determined by (|]), and one 
can use this information to construct the space-time solution of the Einstein equations produced by the initial data 
as a Taylor expansion in t, up to the appropriate order. One gets, using the usual Regge- Wheeler notation the 
following expansions for the £ = 1 components of the metric, 



5 



h ,__» (£±^E"M EMp» +0( <.) pa, 

A' = 3 (^ + v^W)^v^W pt + 0((3) (M) 
)ii= _3 (^ +v ^Wf pi + 0(( , ) (25) 

all the other components being 0(t 3 ). 

The components of the gauge vector generating the first order gauge transformation that makes the initial t = 1 
data purely 0(t 3 ) are, 

£ t =M cos6 (26) 
Cr=Micos9 (27) 
^ = -M 2 sm9 (28) 
6 = 0, (29) 

where, 

P (-3 Vr - 2 MrM - 5 r 3 / 2 M + 2 M 2 Jr + 2 r 5 / 2 + 2 \/r-2Mr 2 ) 
• M ° = 4VrM(r-2M) (30) 

P (-6 Vr - 2 MrM + 2 r 5 / 2 - 8 r 3 / 2 M + 8 M 2 ^ + 3 - 2 MM 2 + 2 Vr-2Mr 2 ) t 
Ml = A^M (31) 

P (-3 Vr - 2 MrM - 5 r 3 / 2 M + 2 M 2 J¥ + 2 r 5 / 2 + 2 Vr-2Mr 2 ) t 

Performing a first order gauge transformation with this generator, one eliminates the first order £ = 1 component 
of the metric. Therefore the leading terms in the initial data become second order. The latter have two contributions, 
both of I — 2 multipolar order. One contribution simply comes from the expansion to second order of the initial 
data generated via the conformal approach. The other contribution comes from the fact that the first order gauge 
transformation we just performed has second order pieces of the form, 

' , 'MP?(l)> !y n pv^(l)TtJ- \°°) 
— -^9\xv,a ,X £(1)£(1) — fffiAxT £(l)£(i)>f ^9\v,a £(l)£(i)>p — <7cta£(1) 

The contribution to the second order Zerilli function coming from the nonlinear terms in the gauge transformation 
can be found after a straightforward but tedious calculation. The relevant contributions to the second order I = 2 
metric at t = 0, are, 



(2) R 



2P 2 R 2 



3M 2 r 3 

\2 



(2) _2 P 2 R 2 (R-Mf 
2 3M 2 r 3 (r-2M) 
(2) _ 1 P 2 R 5 / 2 (R-M) 
1 ~ 3M 2 r 3 / 2 (r-2M) 
(2) lP 2 ^ 

3 M 2 r 3 1 J 



where i? = (y^ + \Jr - 2M) 2 /A. 

Another contribution to this initial data comes from the second order corrections to the conformal factor, calculated 
above. These contributions to the second order I = 2 metric all vanish, except for, 

M\ 3 ,<2),^R 2 



(2) P/ 2 = 4P 2 (l + ^) tf\ R y. 
WK = WH 2 , (35) 



G 



and all first order time derivatives vanish. 

The Zerilli function in an arbitrary gauge is given by |l4 



il>(t,r) 



r(r - 2M) 
3(2r + 3M) 

^2 



^H 2 - 



d r~3M 



Or 



r-2M 



( 2 ) K +(r-2M)[^^~^^h 1 
or r z 



2r + 3M 

Substituting the above components yields the following initial data for the Zerilli function 

(3M + 8R)R 2 



ij)(t = 0,r) 



Ar(2M + 3R) , m , „, (4r - %M)R x l 2 r 3 l 2 ( ddtf* (R) 

2 ln„ 01^ 1 oln„ 1 o ^2 (-«) 1 ; 



6M 2 r(2r + 3M) 3(2r + 3M) 



3(2r + 3M) 



9r 



(36) 



(37) 



0t 



= 0. 



t=o 



For the conformal factor we computed in (|l3|), this results in, 



tp(t = 0, r) 



P 2 R b 



128 



M 
56 R 



6113 M 2 32607 M 3 



M(M + 2R) 4 (M 2 + 10MR + AR 2 ) \ 3 

3882 M 4 235 M 5 691 M 6 26679 M 7 
+ 25 i? 4 ~ T6""R5~ ~ T6""R6" ~ 1600 R 7 ~ 800 "R8 
7P 2 M(M + 12R)(M + 2R) 2 (\n(2M) - ln(M + 2i?)) 



25 i? 2 100 R 3 
2047 Af 8 2327 M 9 " 



19200 i? 9 



+ 



40i? 3 (M 2 + 10Mi? + 4i? 2 ) 



(38) 



(39) 



D. Physical validity of the perturbative treatment 

Inspecting the initial data for the Zerilli function one finds that it behaves in a non-conventional way. Here is where 
we note significant differences with the case of single spinning holes §. 

The first unusual thing one notices is that the initial data goes asymptotically to a constant value for R — > oo. This 
is different from the data for the "close limit" of two black holes (momentarily stationary or boosted) where the Zerilli 
function goes to zero at infinity. The root of this problem can be traced down to the falloff conditions of the extrinsic 
curvature. In all other cases in which the "close limit" approximation has been applied, the extrinsic curvature falls 
off as 1/r 3 at infinity. For the case in this paper, it decreases as 1/r 2 (otherwise the momentum would vanish). This, 
in particular, leads to falloff conditions in the gauge vectors we use to eliminate the t = 1 even pieces of the extrinsic 
curvature. In turn, the falloff condition of the gauge vectors influences, via the quadratic terms in the second order 
gauge transformation, the behavior of the Zerilli function we obtain at the end of the process. 

Does the appearance of a Zerilli function that does not vanish at infinity indicate something problematic in itself? 
We do not seem to see any difficulty in evolving the problem in this case. Our evolution code evolves a "slab" 
region consistent of the initial data and its domain of dependence, therefore we do not need to specify any boundary 
conditions. If one observes the behavior of the Zerilli function as a function of time for a fixed (finite) value of the 
radius, it starts having a constant value followed by a quasinormal ringing and a power law tail decrease towards zero. 
That is, the constant behavior at infinity translates itself in a certain behavior at the beginning of the waveform, it 
does not leave any visible effect after the ringdown and power law tail behavior. The radiated energy and all observable 
quantities at infinity (even to second order in perturbation theory, see section III.B of j^]) are not determined by the 
Zerilli function itself but by its time derivative, for which the initial data indeed goes to zero for large values of the 
radius, and therefore no problem in the evaluation of physical quantities is present. 

Another aspect that could cause concern is the nature of the gauge transformation considered. The gauge trans- 
formation is well behaved in any finite value of r* (that is, any point of the exterior of the black hole excluding 
the horizon and spatial infinity). Since in order to compute the radiated energies and waveforms we will never need 
information from either the horizon nor infinity, the gauge transformation is well defined in all relevant points for our 
calculation. 
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III. RESULTS AND CONCLUSIONS 



We have numerically evolved the Zerilli equation with the initial data presented above and computed the energy 
and waveforms for the "relaxation" of a single boosted Bowen-York hole to a boosted Schwarzschild black hole. The 
results for the energy can be summarized by a simple formula which, to leading order in P reads, 

V M ADM J \MadM J 

where the prefactor was computed numerically 0. 

Figure ^ shows the radiated energy as a function of the momentum. We see that for values of the momentum close 
to P/Madm ~ 0.4 the total radiated energy by the "relaxation" of a single boosted Bowen-York black hole to a 
Schwarzschild black hole is similar to the total radiated energy in a close limit collision of Bowen-York black holes 
|jTTf and figure ||. One could therefore be tempted to say that for values of P/Madm > 0.4 one should stop using these 
families of initial data. A puzzling element is that we have already evolved these families of initial data in collision 
situations and compared with actual non-linear integrations of the Einstein equations fl^ , |Tl|| and we know that these 
families of initial data radiate less in collisions than the values we are predicting here for each individual hole, at least 
for close separations. The results for these collisions are recollected in figure |[ 

What is going on? One has to keep these results in perspective, since it is easy to get carried away and believe that 
perturbation theory should work way beyond where it was meant to do so. In the calculation of the present paper we 
find that the radiated energy goes as P 4 . This is good, since our perturbative parameter is P and therefore this means 
that the corrections are small. In the case of colliding black holes the radiated energies contain terms that go as (see 
reference jy]]) L 4 , P 2 L 2 and PL 3 , where L is the separation of the black hole centers in the conformally related flat 
space. If one simply increases the value of P keeping L constant, it is obvious that the contribution we consider in 
this paper will quickly dominate. However, one is clearly pushing things beyond the realm in which these calculations 
were meant to be reliable. One is essentially forcing a higher order term in perturbation theory to a regime in which 
it is larger than the lower order terms! 

Therefore the conclusion of this paper has to be read in the following way: the terms coming from the "relaxation" 
of a conformally flat hole to a usual boosted slice of Schwarzschild are higher order in perturbation theory than the 
ones one obtains in a collision. Because of this fact, they grow fast with the perturbative parameter and perturbation 
theory breaks down early for the estimation of the involved energy in the calculations of the current paper. The 
breakdown occurs earlier for a single hole than for a collision of holes (at least computed in the center-of-momentum 
frame) . 

Connected with the latter point is another interesting insight gained from the analysis of this paper. It has to do 
with the choice of frame used to describe collisions of black holes in perturbation theory. The situation is illustrated in 
figure One would expect these two collisions to be physically equivalent. However, if one considers conformal black 
holes and expands in perturbation theory, the collision at the bottom will contain terms that behave exactly like those 
we discuss in the current paper, whereas the top collision does not. If we power count, for the bottom collision, the 
extrinsic curvature has leading terms that behave like P and are I = 1 and the conformal factor has terms that behave 
like L 2 plus terms that behave like P 2 at leading order. The energy, being quadratic in the fields, will generically 
contain terms that behave like L 4 , P 2 L 2 and P 4 , the latter being the terms we encountered in this paper. However, 
the top collision contains terms L 4 , PL 3 and P 2 L 2 , as discussed in reference jll]] (the extrinsic curvature goes as PL 
and the conformal factor as L 2 at leading order). Therefore if one were to consider "close" black holes and were to 
consider the radiated energy as a function of P, as we have done above, one would encounter that the perturbative 
predictions — at the order considered — will differ significantly as soon as the P 4 terms start to grow. The moral from 
this paper insofar as the choice of the origin is: perturbation theory breaks down quickly as a function of momentum 
for situations with net linear momentum, it is best to analyze collisions set up in the center-of-momentum frame. 



2 The ADM mass of the slice depends on the value of the momentum. For these calculations we are using the zeroth order 
approximation to the ADM mass, which is constant. It is known from full numerical calculations that for P/Madm < 0.6 the 
discrepancy between the zeroth order approximation and the full value is less than a 1% (see figure 1 of reference pot). This can 
also be seen from the perturbative calculation of the mass, which one can obtain from the conformal factor we computed, and 
yields Madm — M + — j f° r instance, for the "puncture" case corresponding to the conformal factor of equation 
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As expected, given the nature of the Zerilli equation, the waveforms that one obtains from the "relaxation" just 
behave like quasinormal ringing. Figure ^ shows the waveform of the decay. The form is that of a typical ringdown, 
and it therefore takes an amount of time of the order of the light-crossing time of the black hole size to decay. 

Summarizing, as in the case of spinning Bowen-York black holes, one has extra radiation present in the initial data, 
that grows with the value of the momentum. In evolutions of binary black hole collisions, one can either wait long 
enough for this energy to be "flushed away" from the system, or restrict oneself to values of the momentum that are 
small enough that the extra energy is small respect to the total energy produced in the collision. The latter is the 
only option in the case of "close limit" collisions. Another conclusion is that perturbation theory breaks down quickly 
as a function of the momentum of the holes for single boosted holes (and collisions of black holes not computed in 
the center-of-momentum frame) and therefore cannot be used to reliably estimate the "energy content" of each hole 
in a regime that might be of interest for the momenta and energies relevant for black hole collisions. 
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FIG. 1. The energy radiated in the "relaxation of a single, boosted Bowen-York black hole to a Schwarzschild black hole, as 
a function of the momentum. At about P/Madm ~ 0.4 the radiation is equivalent to the total energy radiated in a the "close 
limit" collision of two boosted black holes. 
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FIG. 2. The energy radiated in the head on collision of two boosted black holes as calculated in reference 10. The collision is 
for two black holes separated by a Misner parameter fir, = 1.5. Depicted are the results of first and second order perturbation 
theory and the full numerical results of the NCSA-Potsdam-WashU group. We see that the collisions never radiate more than 
1% of the ADM mass in gravitational waves. 
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FIG. 3. Two black hole collisions that should be physically equivalent, but that are significantly different from the point of 
view of the perturbative treatment involved in the "close limit approximation" . 
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FIG. 4. The time derivative of the Zerilli function as a function of time. The square of this quantity is proportional to the 
energy flux, and therefore characterizes the "waveform" of the gravitational radiation from the decay of a Bowen-York boosted 
black hole into Schwarzschild. 
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The conformally flat families of initial data typically used in numerical relativity to represent 
boosted black holes are not those of a boosted slice of the Schwarzschild spacetime. If such data are 
used for each black hole in a collision, the emitted radiation will be partially due to the "relaxation" 
of the individual holes to "boosted Schwarzschild" form. We attempt to compute this radiation by 
treating the geometry for a single boosted conformally flat hole as a perturbation of a Schwarzschild 
black hole, which requires the use of second order perturbation theory. In this we attempt to mimic 
a previous calculation we did for the conformally flat initial data for spinning holes. We find that the 
boosted black hole case presents additional subtleties, and although one can evolve perturbatively 
and compute radiated energies, it is much less clear than in the spinning case how useful for the 
study of collisions are the radiation estimates for the "spurious energy" in each hole. In addition 
to this we draw some lessons on which frame of reference appears as more favorable for computing 
black hole collisions in the close limit approximation. 
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There is significant interest in obtaining waveforms for the gravitational radiation produced in the collision of black 
holes. Progress is being made on this problem both using supercomputers |l|] and perturbative calculations ||. One of 
the open issues is what families of initial data are appropriate to represent the collision of two black holes, especially 
when the latter are not far away from each other. The state of the art of numerical simulations suggests that for 
some time we may not be able to start simulations with the black holes at a sufficiently large separation, such that 
^ ! one can assume a simple linear superposition will work. This leaves open the issue of how much "spurious radiation" 
I , is one introducing in the various proposals for superpositions in the non-linear regime. Bowen and York || (and with 
J~ I ■ a different set of boundary conditions more recently Brandt and Briigmann || ) studied the problem of giving initial 
data for boosted and spinning holes in such a way that a superposition is possible. They assume the spatial metric 
is conformally flat and as a consequence one can superpose the extrinsic curvatures for the two holes and still solve 
' the momentum constraint. One then proceeds to find a conformally flat spatial metric for the superposed holes by 
solving a nonlinear elliptic partial differential equation. The procedure achieves superposition at the price of assuming 
conformal flatness of the three metrics, which is not generically possible, and more importantly, is not possible in 
situations of interest. For instance a single spinning black hole, described by the Kerr solution, is not known to admit 
conformally flat spatial slices (3|Q| . 

A similar situation develops for the case of a single boosted black hole. The initial data constructed by Bowen and 
York or Brandt and Briigmann do not correspond to those one would find on a boosted slice of the Schwarzschild 
spacetime. In this paper we will refer to these families as "conformal boosted black hole" Therefore if one evolves 
these families of data one should find that the black hole "settles down" to a Schwarzschild form through the emission 
of gravitational radiation. The original purpose of this paper was to study the emitted radiation by treating the 
conformal boosted hole as a perturbation of a Schwarzschild black hole. In a previous paper we had carried out a 
similar discussion for the case of a single spinning conformally flat hole ^ . We will see that the boosted case is much 
more subtle than anticipated. We will be able to evolve the spacetime, but questions will remain about the usefulness 
of the results obtained, at least for the original purpose of gaining insight into the spurious radiation content of data 
for black hole collisions for interesting ranges of parameter values. 

The organization of this paper is as follows. In the next section we will discuss perturbations of a boosted black 
hole. We recall that the £ = 1 even modes are pure gauge and therefore all the physics of interest takes place in the 
£ = 2 modes, which can be treated easily up to second order in perturbation theory. In section III we will discuss 
the perturbative evolution and the amount of radiation produced. We end with a discussion of the results and their 
implication for the choice of frame of reference one uses in perturbative evolutions of black hole collisions with linear 
momentum. 
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II. CONFORMAL BOOSTED BLACK HOLE AS A PERTURBATION OF SCHWARZSCHILD 



A. Initial data in the conformal approach 

The families of initial data that we will consider in this paper are obtained via the "conformal approach" to the 
initial value problem in general relativity. In it, one assumes the metric to be conformally flat g a b = 4- 'Sab and defines 
the conformal extrinsic curvature K a b — <f> 2 K a b- In terms of these variables the initial value constraint equations 
(assuming maximal slicing TrK — 0) read, 

V a K ab = (1) 
1 K ab K >, 

where all the derivatives are with respect to flat space. 

One can construct || solutions to the first set of equations (momentum constraint) for a single black hole centered 
at R — 0, with linear momentum P a , 

■~ 3 

K ab = [ 2P (anb) - {Sab - n a n b )P c n c ] . (3) 

where R is a spherical radial coordinate and n a a radial unit vector, and both are defined in the fiducial flat space 
that one obtains setting the conformal factor to unity. 

Without loss of generality we may assume that P a points along the positive z-axis, and has magnitude P. If we 
write K a b in spherical coordinates, the only non vanishing components are 

-» 3P 
K RR = -^j cos(0) 

^ 3P ~ 

K m = -^ sin ( 9 ) = k or 

— 3P 

Kee = — — cos(0) 

K u = sm{6) 2 K ee (4) 

If we write the extrinsic curvature in terms of tensor spherical harmonics @] , we see that it consists of a pure £ = 1 
even term. This is reasonable, since the presence of momentum in an initial slice configuration is determined by the 
presence of a "dipole" term asymptotically, that makes the ADM integral 

p i = ^zl K ijd 2 S j (5) 



nonvanishing. 

From (0), we have, 



8tt 



KabK =^l( 1 + 2cos20 ) ( 6 ) 



We now need to solve (^), and this involves imposing boundary conditions on <f>. One possibility is that given by 
Bowen and York ||[]. In this case one chooses a certain constant a, solves (||) for R> a, and imposes 

ff + ^ = 0for*=a, (7) 

and 



1 Alternatively, we may use a Brill-Lindquist type boundary condition ^ and as a result one obtains the "puncture" solutions 
off). 
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<f> > , lim 4> = 1. (8) 

R— »oc 

Equation (Q) implies that the inner boundary is an extremal surface, as was shown in reference || and leads to a 
well posed elliptic problem when used with equations (||J|) . One can also use the boundary condition (]?]) to generate 
initial data that obeys an isometry condition, provided one chooses an extrinsic curvature — different from the one we 
are choosing here — that satisfies the isometry condition. This was explored in great detail by Cook || in his numerical 
work constructing solutions for the conformal factor for multiple black holes. The symmetrization procedure yields 
non-trivial results even for a single black hole as we are considering here We have chosen, for simplicity, not 

to symmetrize the extrinsic curvature. Experience has shown that symmetrization does not change significantly the 
amount of radiated energy in head-on black hole collisions Jl2| . Therefore our choice should not crucially influence 
the central conclusions we are attempting to obtain about the conformally flat black hole solutions. 

Even for the simple form of (^), in general we cannot solve @) exactly, and one has to resort to numerical methods. 
In the present analysis, however, we will be interested in solutions for "small" P. Since for P = the solution is 

(O) =1 + | (9) 

we may solve (||) iteratively replacing 

= ^(o) (r) + p 2 ^ 2 ) (i?, 9) + p 4 0( 4 ) (r, e) + ■ ■ ■ (io) 



in and imposing 



(i) 1 

+ — = for R = a, lim W = 0. for i ^ 0, (11) 



as 



dR 2R R^oo 
The form of (||) further suggests that we expand (j>^ in Legendre polynomials Pi(cos9), so that (j)^ may be written 

<j) = <j) W (R) + P 2 U> i o ) {R)Po(cos9)+<t> i 2 ) (R)P 2 (cos9)] +0(P 4 ). (12) 
Solving for the coefficients taking into account the boundary condition (Q) we get, 

^(R)- {R + a)5 - R5 ■ 13 



32a 2 (R + a) 5 512aR 
A2) (J} s _ 75i? 6 + 291aR 5 - 650a 2 i? 4 - 3400a 3 i? 3 - 4800a 4 R 2 - 2925a 5 R - 669a 6 
02 ( )_ 400i?3(i? + a)5 

+ ^l£_[i n (2a) - ln(i? + a)} + (13) 
20i? 3L y ' 128i? 3 V ' 

Notice that 4>^ {R) falls off only as R~ 2 for large R. 

For completeness, we also present the solution one obtains if one chooses the "puncture" boundary condition 
considered by Brandt and Briigmann g recently, 

,( 2);n , M 4 + 10 M 3 R + 40 M 2 R 2 + 80 M R 3 + 80 R 4 , AS 

^ {R) = SM(2R + Mf (W) 

l(2 ) /r ,. 120 R 5 + 768 R 4 M + 1078 R 3 M 2 + 658 R 2 M 3 + 189 RM 4 + 21 M 5 , x 

fflo '(R) = f (15) 

20 J? 2 (2R + M) 

(ln(2 R + M)- ln(M)) M 

40 R 3 

but we will only consider the solution (|l^ ) in the rest of this paper. 

To obtain an initial data set appropriate for a perturbative evolution we proceed as follows. First, in an ADM type 
decomposition, we choose our shift functions N a = 0, so that we have, 

ds 2 = g ab dx a dx b - N 2 (dt) 2 

^ = -2NK ab (16) 
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Next, we assume that on our initial slice, at t = 0, the 3-metric g a b and extrinsic curvature K a b are given by 
the above conformal flatness construction. Namely, we have g a b — <t> A 5 a b and K a b = (f)~ 2 K a \,, with K a b given by 
(f|), and <j> given by (|l2|). (This, of course, ensures that our initial data satisfies the constraint equations). We 
next change from the conformal spherical radial coordinates R, to a "Schwarzschild" radial coordinates r, with 
R = (y/r + y/r — 2M) 2 /4, where M = 2a, and choose N = yl — 2M/r. It can be checked that with these choices, for 
P = the extrinsic curvature vanishes and we recover the Schwarzschild metric in the usual Schwarzschild coordinates. 
This last expression for the metric in (r, 9, <fr) coordinates (or, rather, the initial data) has, therefore, the appropriate 
form for a perturbative treatment of the Regge - Wheeler type, as extended to second order in |jic|| . 

The final form for the initial data that results from this construction has thcfollowing multipolar components: to 
zeroth order in the momentum, we only have £ = components; to first order in the linear momentum, only £ = 1 
contributions; to order 0{P 2 ), we have £ = 0, 2 contributions. All contributions are even-parity and are analyzed in 
more detail in the following Sections. 



B. Multipolar decomposition of the initial data : the I = 1 contributions 



The contribution to the three-metric of zeroth order in P is just the Schwarzschild solution. The first apparently 
non-trivial contribution is given at 0{P l ) and corresponds as we discussed in the previous subsection to an £ — 1 
multipole. Let us analyze the £ = 1 perturbations of a spherically symmetric spacetime. In order to do this we use 
the traditional Regge- Wheeler decomposition. One starts with a background metric written as, 



g$dx»dx u = -(1 - 2M/r) dt 2 + (1 - 2M/r)- 1 dr 2 + r 2 d9 2 + r 2 sin 2 Bdip 2 . 
For axisymmetric perturbations, the general £ = 1 even parity terms take the form 



(17) 



hi 

1.(1 



(1 - 2M/r)H (t,r) cos 6 
Hi(t,r) cos 9 

(1 - 2M/r)- 1 H 2 (t,r) cost 

—ha(t, r) sin 9 

—hi(t, r) sm.9 

r 2 K(t, r) cos 9 

r 2 sin 2 9K(t, r) cos 9 



(18) 



The first order metric perturbation coefficients are not uniquely defined, but may be changed by "gauge transfor- 
mations" of the form 



h^v — hp V 3m^!pC(i) 9np£(iyv 9pv£(xyn 

where the gauge 4-vector £^ is arbitrary, except for the requirement of axisymmetry. 
In particular, the £ = 1 even parity coefficients transform as 



(19) 



H (t,r 

H 2 (t,r 

h (t,r 
hi(t,r 
K(t,r 



H (t,r) 

Hi(f,r) 
H 2 (t,r) 
h Q (t,r) - 
hi(t,r) - 



2 dM (t,r) 2M 

dt r(r-2M)~ 
r -2M dM (t,r) r dMi{t,r) 



-Mx^r) 



r 

2M 



r{r - 2M) 
r-2M 



dr 

M x {t,r)-2 



r - 2M dt 
dMx{t,r) 



M (t,r) - r 
Mx(t,r) - r 



dr 

2 dM 2 (t,r) 
dt 

2 dM 2 (t,r) 



r - 2M 

K(t,r) ~ -Mx(t,r)+2M 2 (t,r) 



dr 



(20) 
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where the functions Mi are arbitrary. 

One can use this gauge freedom to go to a restricted gauge in which ho = hi = K = 0. This gauge is not completely 
fixed. One still can perform gauge transformations of the form, 

where f(t) is arbitrary. 

An interesting result, is that in this gauge it is straightforward to find the general solution of the linearized Einstein 
equations for £ = 1 even parity perturbations. The result is, 

1 (r-2A/) 2 dt 

H 1 F | ^ d2p i 

3(r - 2M) 2 1 3M(r - 2M) 2 dt 2 

H ^V^y Fl{tl (22) 

where F\(t) is an arbitrary function. Remarkably, one can show that this solution is pure gauge. Choosing f(t) = 
—Fi(t)/(6M) in (|2l| ) leads to vanishing gauge transformed quantities. 

We therefore see that the £ = 1 perturbations are pure gauge. This result was first noticed by Zerilli This was 
to be expected in physical grounds since one could always imagine setting coordinates boosted in such a way that the 
black hole would not move. It has the further implication that we can compute the radiated energy by computing 
the £ = 2 perturbations in a gauge in which the £ = 1 perturbations vanish, and studying their evolution. 



C. Multipolar decomposition of the initial data : the I = 2 contributions 

The relevant £ = 2 perturbations are of second order in our perturbation parameter, P. In principle, when one 
works out second order perturbations of a given metric, the evolution equations one gets have the general form of a 
linear operator (similar to the one that acts in first order giving rise to the Zerilli equation) acting on second order 
quantities, equal to a quadratic "source" term formed with the first order perturbations E3 ]. The source term is 
complicated and is delicate to handle numerically when evolving the perturbations. A place where this was explicitly 
done was for instance in the evolution of boosted black hole collisions to second order p^| . The calculations are 
lengthy and complicated. Fortunately, in our case one can proceed in a different way. We have just shown that 
there is a gauge in which the first order perturbations (the £ = 1 ones) vanish. Therefore in that gauge one can 
write a second order Zerilli equation that is source-free. Moreover, the linear portion of the second order perturbative 
equation is exactly the same as the first order perturbative equation (the Zerilli equation), and this equation can be 
written in terms of quantities that are gauge invariant. We notice, however, that eliminating the first order £ = 1 
terms through a first order gauge transformation introduces second-order changes in the metric that are not a pure 
second order gauge transformation, and must be taken into account. One must then be careful in making sure that 
the initial data for the £ = 2 perturbations used in the Zerilli equation corresponds precisely to that gauge. In other 
words, we need to carry out a first order gauge transformation on the initial data that provides a new initial data 
corresponding to a gauge where the £ = 1 perturbations vanish. Since all perturbations satisfy equations that are 
second order in time, this requires that the £ = 1 terms of the metric, and their first time derivative vanish on the 
initial slice. If we consider (|2(]), this requires that the gauge vector components Mi are such that both the left hand 
sides in (po|), and their first time derivatives, vanish when evaluated at t = 0. This, in principle requires only the 
knowledge of the metric functions on the right hand side of (pp|), and their first time derivative, at the fiducial time 
t — 0. It turns out, however, that we also need second order time derivatives, 

to implement the second order gauge transformation required to obtain the £ = 2 initial data. These second 
order time derivatives may be straighforwardly evaluated from the corresponding Einstein equations for the £ = 1 
perturbations. 

For the particular case in question, the £ = 1 initial data (and first time derivative) is determined by (|]), and one 
can use this information to construct the space-time solution of the Einstein equations produced by the initial data 
as a Taylor expansion in t, up to the appropriate order. One gets, using the usual Regge- Wheeler notation the 
following expansions for the £ = 1 components of the metric, 
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h ,__» (£±^E"M EMp» +0( <.) pa, 

A' = 3 (^ + v^W)^v^W pt + 0((3) (M) 
)ii= _3 (^ +v ^Wf pi + 0(( , ) (25) 

all the other components being 0(t 3 ). 

The components of the gauge vector generating the first order gauge transformation that makes the initial t = 1 
data purely 0(t 3 ) are, 

£ t =M cos6 (26) 
Cr=Micos9 (27) 
^ = -M 2 sm9 (28) 
6 = 0, (29) 

where, 

P (-3 Vr - 2 MrM - 5 r 3 / 2 M + 2 M 2 Jr + 2 r 5 / 2 + 2 \/r-2Mr 2 ) 
• M ° = 4VrM(r-2M) (30) 

P (-6 Vr - 2 MrM + 2 r 5 / 2 - 8 r 3 / 2 M + 8 M 2 ^ + 3 - 2 MM 2 + 2 Vr-2Mr 2 ) t 
Ml = A^M (31) 

P (-3 Vr - 2 MrM - 5 r 3 / 2 M + 2 M 2 J¥ + 2 r 5 / 2 + 2 Vr-2Mr 2 ) t 

Performing a first order gauge transformation with this generator, one eliminates the first order £ = 1 component 
of the metric. Therefore the leading terms in the initial data become second order. The latter have two contributions, 
both of I — 2 multipolar order. One contribution simply comes from the expansion to second order of the initial 
data generated via the conformal approach. The other contribution comes from the fact that the first order gauge 
transformation we just performed has second order pieces of the form, 

' , 'MP?(l)> !y n pv^(l)TtJ- \°°) 
— -^9\xv,a ,X £(1)£(1) — fffiAxT £(l)£(i)>f ^9\v,a £(l)£(i)>p — <7cta£(1) 

The contribution to the second order Zerilli function coming from the nonlinear terms in the gauge transformation 
can be found after a straightforward but tedious calculation. The relevant contributions to the second order I = 2 
metric at t = 0, are, 



(2) R 



2P 2 R 2 



3M 2 r 3 

\2 



(2) _2 P 2 R 2 (R-Mf 
2 3M 2 r 3 (r-2M) 
(2) _ 1 P 2 R 5 / 2 (R-M) 
1 ~ 3M 2 r 3 / 2 (r-2M) 
(2) lP 2 ^ 

3 M 2 r 3 1 J 



where i? = (y^ + \Jr - 2M) 2 /A. 

Another contribution to this initial data comes from the second order corrections to the conformal factor, calculated 
above. These contributions to the second order I = 2 metric all vanish, except for, 

M\ 3 ,<2),^R 2 



(2) P/ 2 = 4P 2 (l + ^) tf\ R y. 
WK = WH 2 , (35) 



G 



and all first order time derivatives vanish. 

The Zerilli function in an arbitrary gauge is given by |lf 



r(r - 2M) 
3(2r + 3M) 
r 2 

+ 2r + 3M 



(2) 



H 2 



d r - 3M l9 , 
r — t —- 1 >K 



dr 



(2) 



K+(r- 2M) 



r-2M 

d 

dr 



Substituting the above components yields the following initial data for the Zerilli function 

(3M + 8R)R 2 



tp{t = 0, r) = P z 



6M 2 r{2r + 3M) 3(2r + 3M) 



r) 



4r(2M + 3R) , (2 ) , m (4r - 8M )R 1 ' 2 r 3 / 2 ( d^ (R) * 

2 (if) - 



3(2r + 3M) 



(9 c 



dt 



0. 



For the conformal factor we computed in (|l3|), this results in, 

P 2 R 6 /128 
M(M + 2R) A (M 2 + WMR + 4R 2 ) \~Y 
3882 M 4 235 M 5 691 M 6 26679 M 7 



V(t = 0,r) = 



M 
56 R 



6113 M 2 32607 M 3 

+ 



25 R 2 100 R 3 
2047 M 8 2327 M 9 



25 i? 4 16 i? 5 16 R 6 1600 i? 7 800 i? 8 
7P 2 M(M + 12R)(M + 2R) 2 (\n(2M) - ln(M + 2i?)) 
40i? 3 (A/ 2 + 1QMR + AR 2 ) ' 



19200 R 9 



(36) 

(37) 
(38) 



(39) 



D. Physical validity of the perturbative treatment 

Inspecting the initial data for the Zerilli function one finds that it behaves in a non-conventional way. Here is where 
we note significant differences with the case of single spinning holes || . 

The first unusual thing one notices is that the initial data goes asymptotically to a constant value for R — > oo. This 
is different from the data for the "close limit" of two black holes (momentarily stationary or boosted) where the Zerilli 
function goes to zero at infinity. The root of this problem can be traced down to the falloff conditions of the extrinsic 
curvature. In all other cases in which the "close limit" approximation has been applied, the extrinsic curvature falls 
off as 1/r 3 at infinity. For the case in this paper, it decreases as 1/r 2 (otherwise the momentum would vanish). This, 
in particular, leads to falloff conditions in the gauge vectors we use to eliminate the 1=1 even pieces of the extrinsic 
curvature. In turn, the falloff condition of the gauge vectors influences, via the quadratic terms in the second order 
gauge transformation, the behavior of the Zerilli function we obtain at the end of the process. 

Does the appearance of a Zerilli function that does not vanish at infinity indicate something problematic in itself? 
We do not seem to see any difficulty in evolving the problem in this case. Our evolution code evolves a "slab" 
region consistent of the initial data and its domain of dependence, therefore we do not need to specify any boundary 
conditions. If one observes the behavior of the Zerilli function as a function of time for a fixed (finite) value of the 
radius, it starts having a constant value followed by a quasinormal ringing and a power law tail decrease towards zero. 
That is, the constant behavior at infinity translates itself in a certain behavior at the beginning of the waveform, it 
does not leave any visible effect after the ringdown and power law tail behavior. The radiated energy and all observable 
quantities at infinity (even to second order in perturbation theory, see section III.B of |l0|]) are not determined by the 
Zerilli function itself but by its time derivative, for which the initial data indeed goes to zero for large values of the 
radius, and therefore no problem in the evaluation of physical quantities is present. This might be at first surprising, 
but in reality the Zerilli function deos noe play any physical role, it is its time derivative what does, and the latter 
does not have any unusual asymptotic behavior. To put it in different terms: one could carry out perturbation theory 
and work out all relevant physics entirely in terms of the time derivative of the Zerilli function (which also satisfies 
the Zerilli equation) and there one would not see any unusual asymptotic behavior. 

Another aspect that could cause concern is the nature of the gauge transformation considered. The gauge trans- 
formation is well behaved in any finite value of r* (that is, any point of the exterior of the black hole excluding 
the horizon and spatial infinity). Since in order to compute the radiated energies and waveforms we will never need 
information from either the horizon nor infinity, the gauge transformation is well defined in all relevant points for our 
calculation. 
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III. RESULTS AND CONCLUSIONS 



We have numerically evolved the Zerilli equation with the initial data presented above and computed the energy 
and waveforms for the "relaxation" of a single boosted Bowen-York hole to a boosted Schwarzschild black hole. The 
results for the energy can be summarized by a simple formula which, to leading order in P reads, 

V M ADM J \MadM J 

where the prefactor was computed numerically 0. 

Figure ^ shows the radiated energy as a function of the momentum. We see that for values of the momentum close 
to P/Madm ~ 0.4 the total radiated energy by the "relaxation" of a single boosted Bowen-York black hole to a 
Schwarzschild black hole is similar to the total radiated energy in a close limit collision of Bowen-York black holes 
|L2| and figure ||. One could therefore be tempted to say that for values of P/Madm > 0.4 one should stop using these 
families of initial data. A puzzling element is that we have already evolved these families of initial data in collision 
situations and compared with actual non-linear integrations of the Einstein equations fl6| , |l2|| and we know that these 
families of initial data radiate less in collisions than the values we are predicting here for each individual hole, at least 
for close separations. The results for these collisions are recollected in figure |[ 

What is going on? One has to keep these results in perspective, since it is easy to get carried away and believe that 
perturbation theory should work way beyond where it was meant to do so. In the calculation of the present paper we 
find that the radiated energy goes as P 4 . This is good, since our perturbative parameter is P and therefore this means 
that the corrections are small. In the case of colliding black holes the radiated energies contain terms that go as (see 
reference |fL2[ ) L 4 , P 2 L 2 and PL 3 , where L is the separation of the black hole centers in the conformally related flat 
space. If one simply increases the value of P keeping L constant, it is obvious that the contribution we consider in 
this paper will quickly dominate. However, one is clearly pushing things beyond the realm in which these calculations 
were meant to be reliable. One is essentially forcing a higher order term in perturbation theory to a regime in which 
it is larger than the lower order terms! 

Therefore the conclusion of this paper has to be read in the following way: the terms coming from the "relaxation" 
of a conformally flat hole to a usual boosted slice of Schwarzschild are higher order in perturbation theory than the 
ones one obtains in a collision. Because of this fact, they grow fast with the perturbative parameter and perturbation 
theory breaks down early for the estimation of the involved energy in the calculations of the current paper. The 
breakdown occurs earlier for a single hole than for a collision of holes (at least computed in the center-of-momentum 
frame) . 

Connected with the latter point is another interesting insight gained from the analysis of this paper. It has to do 
with the choice of frame used to describe collisions of black holes in perturbation theory. The situation is illustrated in 
figure One would expect these two collisions to be physically equivalent. However, if one considers conformal black 
holes and expands in perturbation theory, the collision at the bottom will contain terms that behave exactly like those 
we discuss in the current paper, whereas the top collision does not. If we power count, for the bottom collision, the 
extrinsic curvature has leading terms that behave like P and are I = 1 and the conformal factor has terms that behave 
like L 2 plus terms that behave like P 2 at leading order. The energy, being quadratic in the fields, will generically 
contain terms that behave like L 4 , P 2 L 2 and P 4 , the latter being the terms we encountered in this paper. However, 
the top collision contains terms L 4 , PL 3 and P 2 L 2 , as discussed in reference JlJ] (the extrinsic curvature goes as PL 
and the conformal factor as L 2 at leading order). Therefore if one were to consider "close" black holes and were to 
consider the radiated energy as a function of P, as we have done above, one would encounter that the perturbative 
predictions — at the order considered — will differ significantly as soon as the P 4 terms start to grow. The moral from 
this paper insofar as the choice of the origin is: perturbation theory breaks down quickly as a function of momentum 
for situations with net linear momentum, it is best to analyze collisions set up in the center-of-momentum frame. 



2 The ADM mass of the slice depends on the value of the momentum. For these calculations we are using the zeroth order 
approximation to the ADM mass, which is constant. It is known from full numerical calculations that for P/Madm < 0.6 the 
discrepancy between the zeroth order approximation and the full value is less than a 1% (see figure 1 of reference |ll[). This can 
also be seen from the perturbative calculation of the mass, which one can obtain from the conformal factor we computed, and 
yields Madm — M + — j f° r instance, for the "puncture" case corresponding to the conformal factor of equation 
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As expected, given the nature of the Zerilli equation, the waveforms that one obtains from the "relaxation" just 
behave like quasinormal ringing. Figure ^ shows the waveform of the decay. The form is that of a typical ringdown, 
and it therefore takes an amount of time of the order of the light-crossing time of the black hole size to decay. 

Summarizing, as in the case of spinning Bowen-York black holes, one has extra radiation present in the initial data, 
that grows with the value of the momentum. In evolutions of binary black hole collisions, one can either wait long 
enough for this energy to be "flushed away" from the system, or restrict oneself to values of the momentum that are 
small enough that the extra energy is small respect to the total energy produced in the collision. The latter is the 
only option in the case of "close limit" collisions. Another conclusion is that perturbation theory breaks down quickly 
as a function of the momentum of the holes for single boosted holes (and collisions of black holes not computed in 
the center-of-momentum frame) and therefore cannot be used to reliably estimate the "energy content" of each hole 
in a regime that might be of interest for the momenta and energies relevant for black hole collisions. 
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FIG. 1. The energy radiated in the "relaxation of a single, boosted Bowen-York black hole to a Schwarzschild black hole, as 
a function of the momentum. At about P/Madm ~ 0.4 the radiation is equivalent to the total energy radiated in a the "close 
limit" collision of two boosted black holes. 
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FIG. 2. The energy radiated in the head on collision of two boosted black holes as calculated in reference 10. The collision is 
for two black holes separated by a Misner parameter fir, = 1.5. Depicted are the results of first and second order perturbation 
theory and the full numerical results of the NCSA-Potsdam-WashU group. We see that the collisions never radiate more than 
1% of the ADM mass in gravitational waves. 
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FIG. 3. Two black hole collisions that should be physically equivalent, but that are significantly different from the point of 
view of the perturbative treatment involved in the "close limit approximation" . 
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FIG. 4. The time derivative of the Zerilli function as a function of time. The square of this quantity is proportional to the 
energy flux, and therefore characterizes the "waveform" of the gravitational radiation from the decay of a Bowen-York boosted 
black hole into Schwarzschild. 
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